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The identification of phenomena able to pin¬ 
point quantum interference is attracting large in¬ 
terest. Indeed, a generalization of the Hong-Ou- 
Mandel effect valid for any number of photons 
and optical modes would represent an important 
leap ahead both from a fundamental perspective 
and for practical applications, such as certifica¬ 
tion of photonic quantum devices, whose com¬ 
putational speedup is expected to depend criti¬ 
cally on multiparticle interference. Quantum dis¬ 
tinctive features have been predicted for many 
particles injected into multimode interferometers 
implementing the Fourier transformation in the 
Fock space. In this work we develop a scalable 
approach for the implementation of quantum fast 
Fourier transform using 3-D photonic integrated 
interferometers, fabricated via femtosecond laser 
writing technique. We observe the quantum sup¬ 
pression of a large number of output states with 
4- and 8-mode optical circuits: the experimen¬ 
tal results demonstrate genuine quantum interfer¬ 
ence between the injected photons, thus offering 
a powerful tool for diagnostic of photonic plat¬ 
forms. 


I. INTRODUCTION 

The amplitude interference between wavefunctions cor¬ 
responding to indistinguishable particles lies at the very 
heart of quantum mechanics. Right after the introduc¬ 
tion of laser amplification, the availability of strong co¬ 
herent pulses allowed to test interference between differ¬ 
ent light pulses^’^, while generation of pairs of identi¬ 
cal photons through parametric fluorescence^ led subse¬ 
quently to the milestone experiment of Hong, On and 
Mandel^“^. Later on, photonic platforms have been 
demonstrated to be in principle capable to perform uni¬ 
versal quantum computing^. 

Recently, multi-particle interference effects of many 
photons in large interferometers are attracting a strong 
interest, as they should be able to show unprecedented 
evidences of the superior quantum computational power 
compared to that of classical devices^^ The main ex¬ 
ample is given by the Boson Sampling^^ computational 


problem, which consists in sampling from the probability 
distribution given by the permanents of the nxn subma¬ 
trices of a given Haar-random unitary. The problem is 
computationally hard (in n) for a classical computer since 
calculating the permanent of a complex-valued matrix 
is a ^P-hard problem, while it can be efficiently solved 
by n indistinguishable photons evolving through an opti¬ 
cal interferometer implementing the unitary transforma¬ 
tion in the Fock space, and detecting output states with 
an array of single-photon detectors. The chance to pro¬ 
vide evidences of a post-classical computation with this 
relatively simple setup has triggered a large experimen¬ 
tal effort, leading to small-scale implementations^^ as 
well as theoretical analyses on the effects of experimen¬ 
tal imperfections^^’^^ and on possible implementations 
including alternative schemes^^’^^. 

In the context of searching for experimental evidences 
against the extended Church-Turing thesis, a Boson 
Sampling experiment poses a problem of certification 
of the result’s correctness in the computationally-hard 
regime^^. The very complexity of the Boson Sampling 
computational problem precludes the use of a brute-force 
approach, i.e. calculating the expected probability dis¬ 
tribution at the output and comparing it with the col¬ 
lected data. Efficient statistical techniques able to rule 
out trivial alternative distributions have been proposed^^ 
and tested^^’^^, but the need for more stringent tests able 
to rule out less trivial distributions has led, and continues 
to encourage, additional research efforts in this direction. 

In particular, an efficient test able to confirm true n- 
photon interference in a multimode device has been re¬ 
cently proposed^^. The protocol relies on the suppres¬ 
sion of specific output configurations^^ in an interferom¬ 
eter implementing an n^-dimensional Quantum Fourier 
Transform (QFT) matrix, withp being any integer. QFT 
is the Discrete Fourier Transform (DFT) applied to the 
vector of amplitudes of a quantum state, and hence its 
matrix representation is also called Fourier matrix. This 
suppression effect is due to granular^^ many-particle in¬ 
terference and is thus able to rule out alternative mod¬ 
els requiring only coarse-grained effects like the ones 
present in Bose-Einstein condensates^^”^^. Indeed, the 
implications of this effect go well beyond the certifica¬ 
tion of Boson Sampling devices. As a generalization of 
the 2-photon/2-modes Hong-Ou-Mandel (HOM) effect. 
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the quantum suppression law is important at a funda¬ 
mental level, while at the practical level it could be 
used as a diagnostic tool for a wide range of photonic 
platforms^'^’^^’^^’^^. 

In this article, we report the experimental observa¬ 
tion of the recently theoretically proposed^^ suppression 
law for QFT matrices, and its use to validate quantum 
many-body interference against alternative non-trivial 
hypotheses resulting in similar output probability dis¬ 
tributions. QFTs have been implemented with an effi¬ 
cient and reliable approach by exploiting the quantum 
version of the Fast Fourier Transform, an algorithm op¬ 
timized to compute the DFT of a sequence. Here we 
implement Quantum Fast Fourier Transforms (QFFTs) 
on photonic integrated interferometers, using a theoret¬ 
ical proposal developed by Barak and Ben-Aryeh^^, by 
exploiting the 3-D capabilities of the adopted integrated 
circuits. The optical circuits are realized by femtosec¬ 
ond laser writing^^’^^, which makes it possible to fabri¬ 
cate waveguides arranged in three-dimensional structures 
with arbitrary layouts^^”^^. The observations have been 
carried out with two-photon Fock states injected into 4- 
mode and 8-mode QFFTs. The peculiar behaviour of 
Fock states compared to other kinds of states is investi¬ 
gated, showing in principle the validity of the certification 
protocol for the identification of real granular n-particle 
interference, which is the source of a rich landscape of 
quantum effects such as the computational complexity of 
Boson Sampling. 


A. Quantum suppression law in Fourier transform 
matrices 


As a generalization of the HOM effect, it has been 
pointed out that quantum interference effects in mul¬ 
timode interferometers may determine suppression of a 
large fraction of the output configurations^^’^^’^^, de¬ 
pending on the specific unitary transformation being im¬ 
plemented and on the symmetry of the input state. In 
particular^^’^^, let us consider a cyclic input, i.e. an 
n—photon Fock state over m = nP modes (for some in¬ 
teger p) where the occupied modes are determined by 
the rule = 5 + (r — 1) with r = 1,... ,n and 

s = The index s takes into account the 

fact that there are possible n-photon arrangements 
with periodicity which simply differ by a transla¬ 

tion of the occupational mode labels. For example, for 
n = 2 and m = 4 there are 2^ = 2 possible cyclic states, 
(1,0,1,0) and (0,1,0,1), while for n = 2 and m = 8 there 
are 2^ = 4 possible (collision-free) cyclic inputs, i.e. the 
states (1,0,0,0,1,0,0,0), (0,1,0,0,0,1,0,0), (0,0,1,0,0,0,1,0) 
and (0,0,0,1,0,0,0,1). 

We consider the evolution of such states through an in¬ 
terferometer implementing the QFT over the Fock space, 
described by the unitary matrix 


^QFT _ 


1 


^/m 


A 2 ^ 

e* ^ . 


(1) 


Such evolution results in the suppression of all output 
configurations not fulfilling the equation 


mod (E A:z,nj=0, (2) 

where ki is the output mode of the photon. An in¬ 
teresting application of suppression laws is to certify the 
presence of true many-body granular interference during 
the evolution in the interferometer, ruling out alternative 
hypotheses which would result in similar output proba¬ 
bility distributions. In particular, in the case of QFT 
matrices, the observation of the suppression law (2) al¬ 
lows to certify that the sampled output distribution is 
not produced by neither distinguishable particles nor a 
Mean Field state (MF)^^. The latter is defined as a single 
particle state \ 2 p^) of the form 


= (3) 

^ r=l 

with a random set of phases Or for each state, being 
\j^) a single-particle state occupying input mode j^. This 
state reproduces macroscopic interference effects, such as 
bunching or boson clouding^^, which makes it an optimal 
test bed for the certification protocol based on the quan¬ 
tum suppression law (see Fig.l). 

It is possible to quantify the degree of violation V = 
^forbidden/A^events of the suppression law as the number 
of observed events in forbidden output states divided by 
the total number of events^^. If a Fock state was injected 
in an interferometer implementing the QFT with no ex¬ 
perimental imperfections, a violation P = 0 would be 
observed. In the case of distinguishable photons there is 
no suppression law, and the violation will be simply the 
fraction of suppressed outputs, each one weighted with 
the number of possible arrangements of the n distinguish¬ 
able particles in that output combination. In the case of 
two-photon states, the weighting factor is 2 for collision- 
free outputs and 1 otherwise, and a degree of violation of 
1/2 is expected (see Sect. IID). On the contrary, in the 
case of two-photon MF, bunching effects occur leading 
to an expected degree of violation of half the weighted 
fraction of suppressed outputs (1/4 for two-photon MF). 
It has been shown that the fraction of forbidden outputs 
is always large^^. Hence, a comparison of the observed 
value of V with the expected one represents an efficient 
way, in terms of necessary experimental runs, to discrimi¬ 
nate between Fock states, distinguishable particles states 
and MFs. 


B. Realization of Quantum Fast Fourier Transform 

Let us now introduce the experimental implementa¬ 
tion of the QFFT. The general method to realize an ar- 
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FIG. 1: Quantum suppression law for Fock states in a quantum Fourier transform. Conceptual scheme of the 
protocol: the possible configurations of n photons at the output of an m-mode interferometer can be divided into two categories, 
unforbidden and forbidden, depending on whether they satisfy or not the suppression condition (2), respectively. The pie charts 
show the expected output statistics with different cyclic input states, i.e. composed by different classes of particles, where green 
and red areas represent events with unforbidden and forbidden outputs, respectively. The injection of a cyclic Fock state in an 
m-mode QFT interferometer results in total suppression of forbidden output states. Cyclic states with distinguishable particles 
show no suppression, so that each output combination is equally likely to occur. A Mean Field state, which reproduces some 
of the features of bosonic statistics by single particle interference and phase averaging, shows suppression with highly reduced 
contrast. 



directional coupler phase directional coupler phase 


FIG. 2: Schematic of the structure of the integrated devices. Internal structure of the 4-mode (a) and 8-mode (b) 
integrated interferometers implementing the QFFT in the Fock space. In the 8-mode case, the Barak and Ben-Aryeh algorithm 
requires an additional relabeling of the output modes (not shown in the figure), namely 2 5 and 4 7, to obtain the effective 

Fourier transformation. The mode arrangement has been chosen in a way to minimize bending losses. The insets show the 
actual disposition of the waveguides in the cross-section of the devices. The modes coupled together in each step {Li) of the 
interferometer are joined by segments. The implemented phase shifts in each step {Pi) are also indicated. 


bitrary unitary transformation using linear optics was in¬ 
troduced by Reck et al.^^, who provided a decomposition 
of a unitary of dimension m as a sequence of m(m — l)/2 


beam splitters and phase shifters. However, in the spe¬ 
cial case of DFT matrices a more efficient method has 
been proposed^^’^^, which takes advantage of the DFT 
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FIG. 3: Experimental apparatus for input state preparation, a) The photon source (F: filter, HWP: half-wave plate, 
L: lens, PBS: polarizing beamsplitter, PC: polarization compensator, SHG: second harmonic generation, PDC: parametric 
downconversion). b) Photon injection (extraction) before (after) the evolution through the interferometer (FA: fiber array). 


symmetries to significantly reduce the number of linear 
optical elements required. Based on the classical algo¬ 
rithm of Cooley and Tukey^^, who first introduced the 
Fast Fourier Transform (FFT) as a more efficient way 
to calculate the DFT, Barak and Ben-Aryeh developed 
a quantum analogue in the linear optics domain, lead¬ 
ing to the concept of Quantum Fast Fourier Transform 
(QFFT). This approach, valid for 2^-dimensional Fourier 
matrices, requires only (m/2) log m beam splitters and 
phase shifters, to be compared with the O(m^) elements 
needed for the more general Reck decomposition, thus en¬ 
hancing the compactness and scalability of the platform 
for a more reliable experimental realization. The overall 
linear transformation on the optical modes implemented 
by the QFFT circuit is naturally equivalent to that of 
the QFT, hence 

II. RESULTS 

A. Realization of 3D QFFT interferometers by 
femtosecond laser writing 

Here we introduce a new methodology for an integrated 
implementation of the QFFT which exploits the 3-D ca¬ 
pabilities of the femtosecond laser-writing technique. 

The sequential structure arising from the decomposi¬ 
tion of the m-dimensional DFT using the Barak and Ben- 
Aryeh algorithm is reproduced by the consecutive layers 
shown in Fig. 2. The complex arrangement of pairwise in¬ 
teractions necessary for the QFFT algorithm cannot be 
easily implemented using a planar architecture. However, 
femtosecond laser writing technique allows to overcome 
this issue exploiting the third dimension, arranging the 


waveguides along the bidimensional sections of the inte¬ 
grated chip. 

The strategy can be outlined as follows (see also Sup¬ 
plementary Note 1): the 2^ modes are ideally placed on 
the vertices of a p-dimensional hypercube; in each step 
of the algorithm the vertices connected by parallel edges 
having one specific direction are made to interact by 
a 2-mode Hadamard transformation, with proper phase 
terms. An optical interferometer implementing this pro¬ 
cedure is thus composed of log 2 m = p sections, each em¬ 
ploying m/2 balanced beam splitters and phase shifters. 

We fabricated waveguide interferometers realizing the 
FFT for m = 4 and m = 8 modes in borosilicate glass 
chips using femtosecond laser micromachining^^’^^. A 
schematic representation of these two interferometers is 
given in Fig. 2. According to the scheme outlined above 
and by exploiting the three-dimensional capabilites of 
the fabrication technique, the waveguides are placed, for 
what concerns the cross section of the device, on the ver¬ 
tices of a 2D projection of the p-dimensional hypercube 
(see also Supplementary S-1). Three-dimensional direc¬ 
tional couplers, with proper interaction length and dis¬ 
tance to achieve a balanced splitting, connect in each 
step the required vertices. The insets of Fig. 2 show, 
at each step which modes are connected by directional 
couplers {Li) and the amount of phase shift that needs 
to be introduced in specific modes (Pi). Phase shifters, 
where needed, are implemented by geometrical deforma¬ 
tion of the connecting S-bends. Fan-in and fan-out sec¬ 
tions at the input and output of the devices allows in¬ 
terfacing with 127-/im spaced single-mode fiber arrays. 
Note that in our device geometry, in each step, the ver¬ 
tices to be connected are all at the same relative dis¬ 
tance: this means that, unless geometric deformations 
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are designed where needed, light travelling in different 
modes does not acquire undesired phase delays. Finally, 
it is worth noting that the geometric construction here 
developed is scalable to an arbitrary number of modes. 


B. Photon generation and measurement 

The two implemented interferometers of m = 4 and 
m = S modes are injected by single-photon and two- 
photon states. Photons are generated by a type II para¬ 
metric down-conversion process occuring in a BBO (beta- 
barium-borate) crystal, pumped by a Ti:Sa pulsed laser 
(Fig. 3). Spectral and spatial selection on the two pho¬ 
tons are performed by means of interferential filters and 
single mode fibers respectively. The indistinguishability 
of the photons is then ensured by a polarization compen¬ 
sation stage, and by propagation through independent 
delay lines before injection within the interferometer via 
a single mode fiber array. Delay lines are used to ad¬ 
just the degree of temporal distinguishability between the 
photons, so as to record HOM interference patterns. Af¬ 
ter the evolution occurring through the interferometer, 
photons are collected by a multimode fiber array. The 
detection system consists of 4(8) single photon avalanche 
photodiodes used for the 4- (8-) modes chip. Further de¬ 
tails on the photon generation and detection scheme are 
provided in the Supplementary Information. 

To test the validity of the quantum suppression law, 
we measured the number of coincidences at each forbid¬ 
den output combination injecting cyclic inputs with two 
indistinguishable photons. We recorded the number of 
coincidences for each output combination as a function 
of the temporal delay between the two injected photons, 
using an electronic data acquisition system able to detect 
coincidences between all pairs of output modes. 


C. One- and two-photon measurements in 
integrated Fourier interferometers 

Due to the adopted configuration for the type II para¬ 
metric down-conversion source, which generates two pho¬ 
tons in two distinct spatial modes, we limited ourselves to 
the collision-free input case, in which two different input 
modes of the chip are injected. The degree of violation 

V of the suppression law could simply be evaluated with 
a counting experiment. Alternatively, the same quantity 

V can be expressed as a function of single-photon input- 
output probabilities and of the HOM visibilities, defined 
as 


= 


Cij Qij 


a 


hj 


where Cij is the number of detected coincidences for dis¬ 
tinguishable photons and Qij for indistinguishable pho¬ 
tons. The subscripts (i, j) are the indexes of the two 


output modes, for a given input state. The degree of 
violation can therefore be expressed as: 


V = 


Aforbidden 

Aevents 

, E 

(bj) forbidden 


-Fforbidden — 

A= E 

(bj) forbidden 




( 4 ) 


where {Pi^j) are the probabilities of having photons 
in the outputs i,j in the case of indistinguishable (dis¬ 
tinguishable) particles. Here P[^j can be obtained from 
single-particle probabilities. 

For the 4-mode device, the full set of ( 2 )^= 36 collision- 
free input-output combinations has been measured by 
recording the number of coincidences at the two outputs 
as a function of the temporal delay between the two in¬ 
put photons. Due to the law given by Eq.(2), we ex¬ 
pect to observe four suppressed outcomes (over six possi¬ 
ble output combinations) for the two cyclic input states 
(1,3) and (2,4). Since distinguishable photons exhibit no 
interference, HOM dips in the coincidence patterns are 
expected for the suppressed output states. Conversely, 
peaks are expected in the non-suppressed output com¬ 
binations. The experimental results are shown in Fig.4- 
a, where the expected pattern of four suppressions and 
two enhancements is reported, with average visibilities 
of Fdips = 0.899 ± 0.001 and Fpeaks = -0.951 ± 0.004 
for suppression and enhancement respectively. For the 
cyclic inputs, we also measured the interference patterns 
for the output contributions where the two photons exit 
from the same mode. These cases correspond to a full 
bunching scenario with n = 2, and a HOM peak with 
V = —1 visibility is expected independently from the in¬ 
put state and from the unitary operation^^. This feature 
has been observed for the tested inputs, where an aver¬ 
age visibility of V = —0.969 ± 0.024 has been obtained 
over all full bunching combinations. 

The existence of a general rule for the prediction of 
suppressed output combinations when injecting a cyclic 
Fock state in a QFT interferometer is due to the intrin¬ 
sic symmetry of the problem, as opposed to the general 
Boson Sampling scenario^^. Suppressed outputs for non- 
cyclic inputs can be predicted by calculating the perma¬ 
nent of the unitary submatrix given by the intersection of 
columns and rows of corresponding to the occupied 

input and output modes respectively. The injection of 
the full set of input states, including the non-cyclic ones, 
has been employed for a complete chip reconstruction. 
The complete set of measured dips and peaks is shown in 
Fig. 4-a, highlighting the symmetry in the Fourier trans¬ 
form pattern. All single-photon and two-photon mea¬ 
surements described have been used for a reconstruction 
algorithm^^ which allows to extrapolate the effective uni¬ 
tary implemented by the interferometer . Single¬ 

photon and two-photon data are exploited to retrieve in¬ 
formation on the moduli of the unitary elements and on 
the unitary phases. The results are shown in Fig. 4-c,d. 
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FIG. 4: Quantum suppression law in a 4-mode FFT integrated chip, a) Complete set of 36 visibilities measured 
for all input-output combinations in the 4-mode chip. For each input-output combination, the measured coincidence pattern 
as a function of the time delay is shown (points: experimental data, lines: best-fit curves). Cyclic inputs, highlighted in the 
picture, exhibit enhancement (green) and suppression (red) on cyclic and non-cyclic outputs, respectively. For all points, error 
bars are due to the poissonian statistics of the events. For each visibility, error bars are obtained through a Monte-Carlo 
simulation by averaging over 3000 simulated data sets, b) HOM visibilities for all 36 input-output configurations. From left to 
right: experimental measured visibilities (Vqfft), visibilities calculated from the reconstructed unitary (Free), and visibilities 
calculated from the theoretical unitary (Vqft). c-d) Real and imaginary part of the reconstructed experimental transformation 

(^QFFT) 


The fidelity between the reconstructed unitary 
and the theoretical Fourier transform was F = 

0.980 ± 0.002, thus confirming the high quality of the 
fabrication process. 

For the 8-mode chip we recorded all the ( 2 ) = 28 two- 
photon coincidence patterns, as a function of the rela¬ 
tive delay between the input photons, for each of the 4 
collision-free cyclic inputs and for one non-cyclic input. 
We used single-photon measurements to extrapolate the 
effective modulus of the implemented unitary, while two- 
photon visibilities allowed us to reconstruct the five fab¬ 
rication phases (specified in Fig.2) using a minimiza¬ 
tion algorithm. The non-cyclic input has been chosen in 
a way to maximize the sensitivity of the measurements 
with respect to the five fabrication phases. The results 
are shown in Fig. 5. The fidelity between the recon¬ 
structed unitary and the ideal 8-mode Fourier 

transform is F = 0.956 ± 0.005. 

D. Observation of the quantum suppression law 

The suppression of events which do not satisfy Eq. (2) 
is fulfilled only when two perfectly indistinguishable pho¬ 
tons are injected in a cyclic input of a perfect Fourier 
transform interferometer. In such a case, we would have 
the suppression of all output states whose sum of the in¬ 


dexes corresponding to the occupied modes is odd. For 
the 4-mode and the 8-mode interferometer, this corre¬ 
sponds to 4 (16) suppressed and 2 (12) non-suppressed 
collision-free outputs (each one given by two possible ar¬ 
rangements of the two distinguishable photons), plus 4 
(8) terms with two photons in the same output, each one 
corresponding to a single possible two-photon path. 

The expected violation for distinguishable particles can 
be obtained from classical considerations. Let us con¬ 
sider the case with n = 2. The two distinguishable pho¬ 
tons evolve independently from each other, and the out¬ 
put distribution is obtained by classically mixing single¬ 
particle probabilities. All collision-free terms are equally 
likely to occur with probability q = 2/m^, while full- 
bunching events occur with probability q' = ql2 = Ijm?. 
The degree of violation Pd can then be obtained by mul¬ 
tiplying the probability q by the number of forbidden 
output combinations. As a result, we expect a viola¬ 
tion degree of Pd = 0-5 for distinguishable two-photon 
states. The evaluation of the expected value for a mean 
field state, which is due to single particle bosonic statistic 
effects, requires different calculations^^. It can be shown 
that for n = 2 the degree of violation is Pmf = 0.25. 

For each of the cyclic input, we have evaluated here the 
violation degree Pobs resulting from collected data. By 
measuring the coincidence pattern as a function of the 
path difference Ax = cAr between the two photons, and 
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FIG. 5: Quantum suppression law in a 8-mode FFT integrated chip, a) Set of 28 measured visibilities, corresponding 
to all collision-free output combinations for the input (2,6) of the 8-mode interferometer. For each output combination, the 
measured coincidence pattern as a function of the time delay is shown (points: experimental data, lines: best-ht curves). Red or 
green backgrounds correspond to dips and peaks, respectively. For all points, error bars are due to the Poissonian statistics of 
the events. For each visibility, error bars are obtained through a Monte-Carlo simulation by averaging over 3000 green simulated 
data sets. Bottom-left inset: average visibilities of dips (red bars) and peaks (green bars) observed for the 4 collision-free cyclic 
inputs [(1,5), (2,6), (3,7), (4,8)]. Darker regions correspond to error bars of ±1 standard deviation, c-d) Real and imaginary 
part of the reconstructed experimental transformation and of the QFT 


thus by tuning their degree of distinguishability, we could 
address the transition from distinguishable to indistin¬ 
guishable particles. The value of Dobs as a function of Ax 
has been obtained as V.. -x PpANf^f /NP-), where 

TV A- and TVg, are the number of measured coincidences 
for a given value of Ax and for distinguishable particles 
respectively. Two different regions can be identified. For 
intermediate values of Ax with respect to the coherence 
length of the photons, the measured data fall below the 
threshold Pd, and hence the hypothesis of distinguish¬ 
able particles can be ruled out. Then, for smaller values 
of the path difference up to Ax 0, true two-photon 
interference can be certified since both hypothesis of dis¬ 
tinguishable particles and mean field state can be ruled 
out. The maximum violation occurring at Ax = 0 delay 
can be evaluated using equation (4). The experimental 
results retrieved from the protocol are shown in the tables 
of Fig. 6, in which we compare the values Pobs(O) with 
the expected values for distinguishable particles Pd and 
for a mean field state Pmf- As shown for our implemen¬ 
tation, the robustness of the protocol is ensured by the 
high number of standard deviation separating the values 
in each comparison, thus unambiguously confirming the 
success of the certification protocol. In conclusion, the 
alternative hypotheses of distinguishable particles and of 


a mean field state can be ruled out for all experiments. 


III. DISCUSSION 

We have reported on the experimental observation of 
suppression law on specific output combinations of a mul¬ 
timode Fourier transformation in the Fock space due to 
quantum interference between photons. The observation 
of the suppression effect allowed to rule out alternative 
hypotheses to the Fock state. The use of a novel imple¬ 
mentation architecture, enabled by the 3-D capabilities 
of femtosecond laser micromachining, extends the scala¬ 
bility of this technique to larger systems with lower ex¬ 
perimental effort with respect to other techniques. At 
the same time, the universality of a generalized HOM 
effect with an arbitrary number of particles and modes 
is expected to make it a pivotal tool in the diagnostic 
and certification of quantum photonic platforms. Boson 
Sampling represents a key example, since the scalability 
of the technique is expected to allow efficient certification 
of devices outperforming their classical counterparts. 

Fourier matrices can find application in different con¬ 
texts. For instance, multiport beam splitters described 
by the Fourier matrix can be employed as building blocks 
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FIG. 6: Measured violations. Observed violations Pobs as a function of the path difference |Ax| = c|At| between the two 
photons. Blue shaded regions in the plots correspond to the cases where the hypothesis of distinguishable particles can be 
ruled out. Red regions correspond to the cases when both the hypotheses of distinguishable particles and mean field state can 
be ruled out, and true two-particle interference is present, a), Data for the 4-mode interferometer. Blue points: input (1,3). 
Red points: input (2,4). Blue solid line: theoretical prediction for input (1,3). Red solid line: theoretical prediction for input 

(2.4) . Black dashed line: theoretical prediction for a QFT matrix, b). Data for the 8-mode interferometer. Blue points: input 

(1.5) . Red points: input (2,6). Green points: input (3,7). Magenta points: input (4,8). Golored solid lines: corresponding 
theoretical predictions for the different inputs. Black dashed line: theoretical prediction for a QFT matrix. Tables: violations 
7^obs(0) at Ax = 0 and discrepancies (in sigmas) with the expected values for distinguishable particles (T>d) and MFs (Pmf), 
for the cyclic inputs of the two interferometers. Pobs(O) are calculated following formula (4), while expected values for the 
other two cases are Pd = 0.5 and Pmf = 0.25. Error bars in all experimental quantities are due to the Poissonian statistics 
of measured events. All theoretical predictions in solid lines are calculated from the reconstructed unitary and a two-photon 
state with indistinguishability a — 0.95. 


for multiarm interferometers, which can be adopted 
for quantum-enhanced single and multiphase estimation 
protocols"^^. This would also allow the measurement of 
phase gradients with precision lower than the shot-noise 
limit^^. Other fields where Fourier matrices are rele¬ 
vant include quantum communication scenarios^^, obser¬ 
vation of two-photon correlations as a function of geo¬ 
metric phase^^, fundamental quantum information the¬ 
ory including mutually unbiased bases^^, as well as en¬ 
tanglement generation^^. 

Methods 

Waveguide device fabrication. Waveguide inter¬ 
ferometers are fabricated in EAGLE2000 (Corning Inc.) 


glass chips. To inscribe the waveguides, laser pulses with 
300 fs duration, 220 nJ energy and 1 MHz repetition 
rate from an YbiKYW cavity dumped oscillator (wave¬ 
length 1030 nm) are focused in the bulk of the glass, us¬ 
ing a 0.6 NA microscope objective. The average depth of 
the waveguides, in the three-dimensional interferometric 
structure, is 170 /rm under the sample surface. The fabri¬ 
cated waveguides yield single mode behaviour at 800 nm 
wavelength, with about 0.5 dB cm“^ propagation losses. 
The central part of the three-dimensional interferometer, 
which includes all the relevant couplers, have a cross sec¬ 
tion of about 50 /imx50 /rm (95 /imx95 /im) for a length 
of 9.0 mm (14.7 mm) in the four-(eight-)modes case. The 
length of each fan-in and fan-out section, needed to bring 
the waveguides at 127 /im relative distance, is 7.8 mm 
(13.2 mm). 
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^ Electronic address: fabio.sciarrino@uniromal.it 
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INTERFEROMETER DESIGN IN THE 2^-MODES CASE 

In the following, we will schematically discuss the procedure to design an interferometer by implementing the Fast 
Fourier Transform (FFT) in the general case of m = 2^ modes, for some integer p. 

• Each optical mode (waveguide) k G [l,m] can be numbered according to its binary representation and can be 
associated to the set 

{bi,b2,...,bp)k 

where 6 ^ = —1 or +1 if the i-th bit of the binary representation of k is 0 or 1 respectively. The ordering of the 
biS is intended to be from the most significant bit bi to the least significant bit bp. 

• The p-element vector ( 61 , 62 , ...^bp) represents the coordinates of the vertices of a p-dimensional hypercube in 

• The optical FFT algorithm^ consists of p steps. The j-th step connects, making them interact, all the couples 
of modes that differ only for the j-th bit. (Proper phase terms are also added and a final relabelling of the 
outputs is needed, but this is not relevant here). 

• The position of each waveguide in the cross-section plane can be defined by suitably projecting in two dimensions 
the vertices of the hypercube. Each step of the EET corresponds to connecting, by directional couplers, the 
couples of modes that have to interact in that step. In particular, such connections corresponds to edges of the 
hypercube with a given direction, all parallel to each other (note that the edges of a p-dimensional hypercube 
are placed along p possible directions). 

• It is worth of noting that with this layout, in each step, modes that have to be connected by couplers are all at 
the same relative distance and the projection of the couplers on the plane are all parallel lines. This means that 
such couplers can be all identical (except for deformations that may be needed for introducing phase terms) and 
waveguides never cross. This avoids parasitic coupling between modes that should stay separated and unwanted 
differences in the optical paths of the mode connections. 


TWO-PHOTON SOURCE 

The generation of two-photon states is performed by pumping a 2 -mm long BBO crystal with a 392.5 nm wavelength 
pulsed laser, with average power of 750 mW. After generation in the crystal, the two photons at 785 nm are spectrally 
filtered by means of 3 nm interferential filters, and coupled into single mode fibers. Polarization compensator and 
spatial delay lines are used before injecting the photons into the integrated interferometers through a single mode fiber 
array. Photons collection and detection is performed with a multimode fiber array and avalanche photodiodes. Finally, 
an electronic acquisition system allows to record two photon coincidences between arbitrary pairs of detectors. Typical 
coincidence rates for each collision-free output combination with distinguishable photons amounted to ^ 70 — 80 Hz 
(for the 4-mode chip) and ^10 — 20 Hz (for the 8 -mode chip). 
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Supplementary Figure S - 1. Planar projections of hypercubes of dimension p up to 4. The numbered vertices correspond to the 
possible position of the waveguides in the cross-section of a 3D interferometer which implements the Fast Fourier Transform. 
Each connection between the vertices corresponds to a directional coupler in a given section of the interferometer. 
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